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We investigate the ground states of antiferromagnetic Mn nanochains on Ni(llO) by spin-polarized 
scanning tunneling microscopy in combination with theory. While the ferrimagnetic linear trimer 
experimentally shows the predicted collinear classical ground state, no magnetic contrast was ob- 
served for dimers and tetramers where non-collinear structures were expected based on ab-initio 
theory. This striking observation can be explained by zero-point energy motion for even numbered 
chains derived within a classical equation of motion leading to non classical ground states. Thus, 
depending on the parity of the chain length, the system shows a classical or a quantum behavior. 



Magnetism is ultimately caused by the spin degree of 
freedom of the electrons. If phase coherence of the elec- 
trons was preserved, the quantum nature of spin would 
potentially allow to encode quantum information in spin- 
tronic devices [1]. Thus, realizing magnetic quantum de- 
vices necessarily involves the understanding of the spin 
of nanoscopic structures on a quantum mechanical level. 
Antiferromagnetic nanostructures are by far not as well 
studied as their ferromagnetic counterparts. This defi- 
ciency lies in the inherent experimental and theoretical 
difficulties which have to be overcome to understand anti- 
ferromagnets. Up to date, even for very simple structures 
such as the one-dimensional antiferromagnetic chain, the 
ground state is unknown. While neutron diffraction of 
one dimensional antiferromagnets often revealed a sim- 
ple, i.e. classical, alternative orientations of the spins [2] 
- called the Neel state - half integer spin chains are for 
example expected to be in a complex entangled ground 
state [3]. The ground state becomes even more complex 
when competing exchange interactions exist, leading to 
magnetic frustration, non-collinear spin structures [4] or 
to correlated ground states predicted by the Anderson 
resonating valence bond model [5] . Geometric frustration 
of the Heisenberg antiferromagnet on a triangular lattice 
is the standard example of a magnetically frustrated sys- 
tem [6-9]. Here we show that antiferromagnetic chains 
display a classical Neel state for odd numbered length 
and an entangled state for even numbered length when 
competing exchange interactions and spin-orbit interac- 
tions are present. Thus, the chains alternate between 
opposite nature of the ground states just by the removal 
or addition of a single atom. 

In the present study, i.e. in atomic Mn chains on 
Ni(llO), the frustration arises from the antiferromagnetic 
coupling within the chain competing with the ferromag- 
netic coupling of the chain atoms to the substrate. Lou- 
nis et al showed theoretically, that this can lead to 
an even-odd effect, where the magnetic structure cru- 
cially depends on the parity of the number of atoms in 



the chain [10]. Sole consideration of the predominant 
antiferromagnetic coupling within the chain leads to an 
antiparallel order of the magnetic moments. Therefore, 
odd- numbered chains exhibit a net magnetic moment, in 
contrast to even- numbered ones. Switching on the weaker 
ferromagnetic coupling between the atoms of the chain 
and the substrate thus acts differently on the two kinds of 
chains. While odd-numbered ones retain their collinear- 
ity and the net moment of the chain aligns with that 
of the substrate, even-numbered chains develop a more 
complex ground state. In a presumed collinear state, the 
total magnetic exchange energy of the Mn chains to the 
ferromagnetic substrate is independent of the direction 
of the Mn moments. It can, however, be lowered when 
a non-collinear spin- structure develops. While magnetic 
exchange energy has to be paid to tilt the spins of the 
even-numbered chain from the ideal collinear state, a net 
spin of the chain develops that points in the direction of 
the substrate magnetization thus giving rise to an energy 
gain due to the exchange with the substrate. 

Mn chains of lengths up to 6 atoms on a Ni(llO) sur- 
face were simulated using the Korringa-Kohn-Rostoker 
Green function method [11, 12] as expressed within den- 
sity functional theory (DFT) taking into account non- 
collinear spin structures [10, 13]. The results were then 
mapped to a classical Heisenberg model in which mag- 
netic exchange energies between first neighbor atoms 
were taken into account [14] . As shown before for chains 
on Ni(OOl) [10], this model catches the important fea- 
tures observed in the ab-initio calculations. 

Figure 1 illustrates the ground states of linear dimer, 
trimer and tetramer chains on Ni(llO) calculated within 
this framework. As expected, the trimer chain is in a 
collinear magnetic state where the net spin of the chain 
aligned to the substrate magnetization. The even num- 
bered chains show a non-collinear magnetic configura- 
tion with moments strongly deviating from the substrate 
magnetization direction. Neglecting spin-orbit interac- 
tion, the magnetic moments of the chains can be rotated 
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FIG. 1: DFT ground states of antiferromagnetically cou- 
pled Mn chains on Ni(llO). While the linear trimer shows 
a coUinear ferrimagnetic order with magnetic moments (blue 
vectors) parallel to the Ni magnetization (red arrow), even- 
numbered chains show a non-collinear magnetic structure. 
Without spin-orbit interaction, the magnetic moments can 
be coherently rotated around the Ni magnetization without 
changing the energy as depicted by the black circles. 

coherently around the direction of magnetization of the 
substrate without affecting the total energy. When spin- 
orbit interaction is taken into account, the rotational de- 
generacy is lifted and the Mn moments prefer an out- 
of-plane orientation, while Ni(llO) is magnetized in the 
surface plane due to shape anisotropy. Thus, two degen- 
erate ground states are predicted. 

So far there has been no experimental report on this 
effect. Currently, only spin-polarized scanning tunnel- 
ing microscopy (Sp-STM) can reveal antiferromagnetism 
on the atomic scale [15, 16]. Low temperature STM has 
been used to investigate the quantum nature of small 
magnetic clusters [17-22]. Furthermore, STM is capable 
of moving adatoms thus offering the possibility of assem- 
bling and probing at the same time [23]. In this work a 
home-built STM operating at 4.2 K and under ultrahigh 
vacuum conditions was used in combination with W-tips 
coated with 10 ML of Fe, 15 ML of Mn, or 30 ML of Co 
for the spin-polarized measurements. 

We first deposited 0.02 ML Mn with the sample held 
at 4.2 K showing primarily single Mn adatoms (see Fig. 2 




FIG. 2: STM images of Mn atoms on Ni(llO). (a) Sample 
with 0.02 ML Mn/Ni(110) deposited at 4.2 K showing mainly 
isolated atoms, (b) Formation of a linear Mn trimer by atomic 
manipulation. 



(a)). Mn chains with the intrinsic nearest-neighbour 
spacing of 2.49 A were then assembled by atomic ma- 
nipulation along the close-packed (110) direction of the 
substrate (see Fig. 2 (b)). The assembly was limited to 
tetramers as longer chains were unstable due to the large 
lattice mismatch between Mn and Ni [24]. Tunneling 
spectroscopy was used to determine the electronic struc- 
ture of the chains revealing no resonances below 1 eV (see 
Supplementary material). Using Sp-STM with in-plane 
magnetized Fe coated tips, we investigated the magnetic 
structure of Mn chains on the atomic scale. For t rimers, 
the measurements revealed a strong spin contrast along 
the chain at a sample bias of 350 mV (see Fig. 3(b)). 
Similarly, the line section along the trimer axis displays 
two minima at the edge adatoms and a maximum in their 
centers. This is in full agreement with the predictions of 
antiferromagnetic odd-numbered chains displaying a sim- 
ple collinear spin structure (see Fig. 1). The dimer (Fig. 
3 (a)) does not show a contrast in agreement with the 
predictions. In the calculated magnetic configuration, 
the projection of the individual Mn magnetic moments 
on the substrate magnetization are identical. Neither, 
the tetramer (see Fig. 3 (c)) displays a strong contrast. 
At a closer look, the line section shows a small depres- 
sion of the d//dU signal at the edge atoms of the chain 
(arrows in Fig. 3 (f)) . This is expected for the tetramer, 
as the edge atoms align more to the Ni substrate mo- 
ments. Thus, all chains show the expected projections of 
the magnetic moment along the substrate magnetization. 




I .CTTT)Mnj 

M— 2nm ► Q Q5 2.0 

lateral position [nm] 

FIG. 3: (a-c) Sp-STS results obtained with Fe-coated tip 
at 350 mV, 6 nA with the same tip on a dimer, trimer and 
tetramer, respectively, (d-f) Linescans along the dashed lines 
in the corresponding images (a-c). 

Dimers and tetramers were suggested to have a non- 
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collinear magnetic structure and as a consequence, they 
should display a more complex spin contrast when using 
a tip magnetized perpendicular to the Ni moments. In 
experiments with different Mn tips of random magnetic 
orientation, Fe tips with arbitrary in-plane orientation 
or Co tips with out-of-plane orientation, we, however, 
never observed such a contrast. This rather suggests that 
the expectation value of the Mn moments perpendicular 
to the magnetization direction of the substrate does not 
vary along the even-numbered chains under the experi- 
mental conditions. Possibly, spin-fluctuations are present 
induced by external perturbations, e.g. temperature or 
inelastic scattering induced by the tunneling electrons. 

We, however, argue in the following that when treat- 
ing the energy of the spin configuration within classical 
equations, we do not predict a classical ground state but 
a non-classical state due to zero-point motion even at 
zero temperature [5] . The complexity of the magnetic in- 
teractions in non-collinear structures with spin-orbit in- 
teractions hinders the use of ab-initio methods beyond 
DFT, e.g. a time-dependent procedure (see for example 
Ref. [25, 26]). Also a model Hamiltonian of a quantum 
Heisenberg system with localized spins and exact diago- 
nalization [17] describing effects like magnetization tun- 
neling [27] cannot be used, as it neglects the itinerant 
nature of the system. Instead, we take a new pathway to 
describe even- numbered chains following ideas of Leggett 
et al. [28] in the framework of the "spin-boson" prob- 
lem. By integrating out the electronic degrees of freedom 
described within DFT, we compute the interatomic ex- 
change interaction and the intraatomic spin-orbit interac- 
tion. The parameters extracted from DFT are then fed 
into classical equations of motion. The resulting equa- 
tions are then treated in the limit of quantum mechanics 
to investigate the spin-dynamics of the system. 

The equation of motion [29] for a magnetic atom with- 
out damping is given by 
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magnetic moment M, and 7 is the gyromagnetic ratio. 
Heff can be determined from the total energy according 
to 
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where Hgff is the effective magnetic field acting on the 



where E is the total energy corresponding to a Heisen- 
berg Hamiltonian which includes magnetic anisotropy 
and in which the magnetic exchange energies between 
first neighbors are considered. For the Mn dimer we find: 

E(e, (j)) = - Ji cos(2(9) - 2J2 cos((9) -h 2Ksin^{0) cos^((/)) 

(3) 

Here, is the polar angle of Mn magnetization to the 
magnetization direction of the Ni substrate and (j) is the 
azimuthal angle of the magnetization direction of the Mn 
atoms with respect to the easy axis. The Ni moments 
were treated as rigid. Ji and J2 are the exchange con- 
stants for Mn-Mn and Ni-Mn exchange, respectively. K 
is the uniaxial anisotropy per adatom of the Mn dimer 
for rotation around the Ni magnetization direction. The 
equation describes the coupled dynamics of and ^ in a 
harmonic potential for small deviations from the ground 
state with 00=0°, 180°. 

From our ab-initio calculations we determine Ji = 
-221.4 meV, J2 = 116.8 meV and K = -0.3 meV which 
result in ^0 = 75° and an easy direction of the Mn mo- 
ments along the surface normal. The energy barrier to 
coherently and adiabatically rotate the moments of the 
Mn atoms is proportional to the anisotropy barrier K 
multiplied by 2sin^(^o) = 0.56 meV. When solving the 
coupled equation of motion in and (j) around the ground 
state, we obtain an eigen-frequency u for the precession 
of the magnetic moments. This precession involves both 
(j) and (see e.g. Ref. [30]). It corresponds to a periodic 
oscillation of the two variables in a local well of the po- 
tential. It can be treated to lowest order as a harmonic 
oscillator. This frequency defines thus a zero-point fluc- 
tuation energy Eq'. 
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A/-4i^cos(20o)(4Ji cos(26>o) + 2J2 cos(6>o) + 4i^cos(26>o) cos2(0o)) (4i^cos(6>o) sin(20o))2 



(4) 



where ^ ?^ 2 is the g- factor and /i^ the Bohr magne- 
ton. Surprisingly, also the exchange interactions enter 
the zero-point fluctuation energy leading to Eo of the or- 
der of 8.9 meV. Thus, in the ground state, the magnetic 
moments have a much larger fluctuation energy than the 



anisotropy barrier and the system can overcome the bar- 
rier. 

Similar equations for the tetramer also reveal, that 
the zero-point motion is large enough to overcome the 
anisotropy barrier (see Supplementary material). In the 
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framework of the "spin-boson" problem, this case rep- 
resents the situation where the matrix element for tun- 
neling between the two localized states (mediated by the 
exchange interaction) is larger than the barrier (medi- 
ated by the spin-orbit interaction). In this pathological 
case of the "spin-boson" problem, the states do not lo- 
calize in one or the other classical state even at T = K 
[28] . It is crucial and very instructive to analyze the pre- 
vious equation: The barrier height and Eo depend on 
both the anisotropy and the exchange energies. While 
the maximal value of the barrier height is limited by the 
anisotropy, the zero-point energy strongly depends on the 
exchange constants. One notices that for typical values of 
the exchange being much larger than the anisotropy, the 
zero-point energy is also much larger than the anisotropy. 
By decreasing Ji and omitting J2 (paramagnetic sub- 
strate), the zero-point energy decreases quickly. Thus at 
the limit of weak interactions between the adatoms and 
a finite barrier due to anisotropy, the system exhibits a 
localized, i.e. Neel, ground state [31]. 

To illustrate the equations of motion, we map the prob- 
lem onto a mechanical analogon. In the limit of infinite 
Mn exchange and absence of spin-orbit interaction all 
antiferromagnetic Mn configurations irrespective of their 
orientation to the substrate magnetization direction are 
degenerate. The mechanical analogon would be two spin- 
ning masses (spinning axis indicated by the blue arrows 
representing the magnetization direction) coupled by a 



(a) 




FIG. 4: Classical analogon of the dimer ground state (a) in 
the limit of an infinite antiferromagnetic exchange of the Mn 
dimers and in the absence of spin-orbit interaction and (b) 
when the antiferromagnetic exchange is finite in the presence 
of a magnetic anisotropy. 



rigid stick (see Fig. 4 a). The exchange interaction to the 
substrate creates an effective field equivalent to gravity 
pulling the two masses down. The ground state is, how- 
ever, continuously degenerate without a barrier to rotate 
the spins. Thus in the quantum limit of the mechanical 
model, the dimer can be in an arbitrary superposition of 
different states leading to a vanishing expectation value 
of the magnetic moment for the individual atom. If we 
assume a finite antiferromagnetic exchange, the dimer 
deviates from a strictly antiparallel orientation in the ex- 
change field of the ferromagnet much alike in the me- 
chanical model, when the two masses are connected with 
springs (see Fig. 4b). When taking into account uni- 
axial magnetic anisotropy, two degenerate ground states 
evolve as indicated in Fig. 4b separated by a shallow po- 
tential barrier, shown in red. These ground states repre- 
sent the degenerate states as determined with DFT with 
spin-orbit interaction. As the vibration of the spring is 
coupled to a rotation of the dimer magnetization due to 
precession, we expect that the mechanical model pop- 
ulates a discrete ground state at T = K, in case the 
barrier between the two ground states cannot be over- 
come by zero-point fiuctuations of the spring [28]. If it 
can, a superposition of the two degenerate ground states 
is the true ground state, as in our case due to the large 
zero-point energy caused by the strong Mn-Mn exchange. 

In conclusion, we have shown that although there is 
a finite magnetic anisotropy energy acting as a poten- 
tial barrier between the degenerate DFT ground states 
in even- numbered chains, a zero-point energy provides a 
mean for fiuctuations between the two degenerate states. 
The latter one is found to be surprisingly large and should 
be considered when describing antiferromagnetic nano- 
objects. Just by adding or removing one atom of the 
chain - changing parity - the system changes its magnetic 
behavior completely and behaves classically. This is due 
to the net spin of the chain coupling to the macroscopic 
magnetization of the Ni substrate. We believe that our 
findings can show a path to create magnetically stable 
antiferromagnetic structures, i.e. to raise the blocking 
temperature. This intriguing result is obtained without 
requiring a quantum Heisenberg model but by treating 
the equation of motion for the magnetic moments as a 
quantum equation. In this respect, we believe that this 
approach is general and could lead to a better under- 
standing of the dynamics of small spin systems. 
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Sample preparation 

All experiments have been performed in ultra high vac- 
uum {p < 10-^° mbar). The Ni(llO) crystal was cleaned 
by repeated cycles of Ar+ ion sputtering (1.5 kV, 30 min) 
and annealing to 970 K until no contaminations were 
found by low energy electron diffraction (LEED) . Atomi- 
cally resolved scanning tunneling microscope (STM) im- 
ages revealed a clean surface with terrace sizes exceeding 
200 nm in width. 

Mn was evaporated in situ from high purity Mn grains 
placed into a Mo-crucible heated by electron bombard- 
ment with typical deposition rates of approximately 
0.1 ML/min. For Mn deposition above room tempera- 
ture the formation of an ordered surface alloy at half 
monolayer coverages of Mn on Ni(llO) has been deduced 
from LEED experiments [1, 2]. This phase is character- 
ized by a c(2 x 2) surface unit cell. Thus, deposition 
below room temperature is indispensable. When keep- 
ing the sample at ~ 170 K during the deposition process, 
no superstructure was observed by LEED. This indicates 
that there is no ordered intermixing between manganese 
and nickel. 

Although for low coverages of Mn [r^ 0.02 ML) one ob- 
serves close-packed monoatomic chains of adatoms over 
a wide range of temperatures, the formation of chains 
was always accompanied by isolated protrusions in the 
nickel substrate as seen by STM. We attribute these pro- 
trusions (apparent height 30 pm) to incorporated man- 
ganese atoms, as their density roughly scaled with the 
amount of deposited Mn and is reduced for lower de- 
position temperatures. This was further confirmed by 
comparing scanning tunneling spectroscopy (STS) mea- 
surements performed on protrusions in the substrate and 
on Mn adatoms exhibiting a similar structure. This, in 
turn, implies that such chains of adatoms are intermixed 
Mn/Ni chains. Therefore, direct deposition of Mn on the 
sample held at cryogenic temperatures (4 K) and atomic 
manipulation is necessary for pure Mn chains. 

As STM tips, chemically etched W- wires were used 
that were cleaned in situ by Ar+ ion sputtering and heat- 
ing to 3700 K. For Sp-STS measurements the cleaned 
W-tips were coated with 10 ML of Fe or 15 ML of Mn or 
30 ML of Co and gently annealed. 

The STM was cooled to 4.2 K using a liquid He bath 
cryostat. Mn chains were assembled by atomic manipu- 



lation, i.e., by positioning the tip on a desired atom, re- 
ducing the resistance of the tunneling gap to 40 Ml and 
moving the tip along the groves of the Ni(llO) surface. 
When returning to the normal tunneling conditions of 
gap resistances larger than 1 MQ, the success of the lat- 
eral manipulation was checked by imaging. Sp-STM was 
carried out with coated tips and modulating the bias volt- 
age by 10 mV at high frequencies (10 kHz range). The 
differential conductance was then recorded as function of 
tip position using a lock-in detection. 



Confined electronic states 

In STS measurements, the differential conductance 
d//dU is proportional both to the local density of states 
(LDOS) of the sample and the tunneling matrix ele- 
ment [3]. By means of a proper normalization of the 
d//dU signal, the LDOS of the sample can be inferred 
[4]. Fig. 1 (a) shows the normalized differential conduc- 
tivities measured at the indicated positions of the cor- 
responding Mn chains and on the Ni substrate. As the 
projected electronic structure of Ni(llO) does not exhibit 
a pseudo gap, the retrieved LDOS is rather constant. 

In contrast, all spectra taken on Mn chains show dis- 
tinct features. For a single Mn atom, there is a broad 
peak at 2.4 V (full width at half maximum ~ 0.8 eV). 
Spectra taken at the center of trimers and tetramers have 
a characteristic substructure separable into several over- 
lapping peaks. Changing the lateral position of acqui- 
sition leads to a change in the relative intensity of the 
sub-peaks but leaves their energy unchanged. Maps of 
the local differential conductance at fixed bias voltage 
show a standing wave pattern (see 1 (b) and (c)). The 
peaks in the d//dU spectra arise from electronic states 
confined to the Mn chains. Remarkably, the onset of the 
lowest peak in energy shifts to lower bias voltages for 
longer chains (from 2.4 V for the single atom to 1.7 V for 
the tetramer). This is in qualitative agreement with a 
simple particle in a box model. 

The results for atomic Mn chains on Ni(llO) are con- 
sistent with similar results previously reported for copper 
chains on Cu(lll) [5] and for gold chains on NiAl(llO) 
[6] . Most importantly, no standing wave pattern has been 
observed at low sample bias (< 1 V) indicating a homo- 
geneous LDOS at these energies. 
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Ab-initio calculations 

The ab-initio calculations were performed using 
density functional theory as implemented in the 
full-potential Korringa-Kohn-Rostoker Green function 
method [7], where the non-collinear treatment was 
used [8] including relativistic effects [9]. The method 
is ideal for investigating nanostructures on substrates in 
a real space embedding approach. A Dyson equation is 
solved where the Green function of interest G is related 
to the Green function of the perfect substrate Go via the 
change of potential V induced by the perturbing nanos- 
tructure: G = Gq -\- GqVG (in a matrix notation). From 
the green function G, the charge and spin densities are 
extracted self-consistently. 

The magnetic exchange interactions Jij mentioned in 
the paper and used in the Heisenberg model as well as in 
the formula for the zero-point energy describe the pair- 
wise exchange interactions within the Mn nano- wires. 
Their values were calculated considering infinitesimal ro- 
tation of the magnetic moments [10]. This methods gives 
a good account of the interactions since the results ob- 
tained with the Heisenberg model are very close to the 
first-principles calculations as shown in Ref. [11]. The 
results are shown in Table I. One notices that the total 
length of the chain has a minor infiuence on the near- 



(b) dl/dU maps at 2.3V 





■ 3 nm ■ 



(c) dl/dU maps trimer 



n ' r 

1 2 3 

sample bias (V) 













j 




i 


r 


1 




2.3( 


)V 



■ 2.5 nm - 



FIG. 1. (a) Normalized d//dU spectra taken on Mn 
nanochains on Ni(llO) acquired at the indicated positions in 
the corresponding topographic image, (b) d//dU maps of a 
dimer, trimer, and tetramer at 2.3 V, showing that this state 
is predominantly located at the edges of the chains, (c) The 
d//dU maps of a trimer taken at 1.82 V and at 2.3 V show 
that the states that correspond to the peaks in the d//dU 
spectra have different spatial distributions. 



est neighbor exchange constants (Jn.n. ~ —200 meV), 
whereas there is a dramatic decrease with increasing dis- 
tance. All exchange constants between second nearest 
neighboring adatoms are positive, which explains that 
these adatoms have moments oriented along the same 
directions. 



TABLE L Exchange constants in units of meV between Mn 
adatoms within different nano- wires. 



Jij 


Dimer 


Trimer 


Tetramer 


Pentamer 


Sextamer 


12 


-221.43 


-207.13 


-179.31 


-187.04 


-191.51 


13 


/ 


+31.50 


-hl9.14 


+ 11.76 


+11.58 


14 


/ 


/ 


-8.77 


-7.17 


-5.05 


15 


/ 


/ 


/ 


+2.4 


+3.26 


16 


/ 


/ 


/ 


/ 


-0.57 


23 


/ 


-207.13 


-191.00 


-172.63 


-180.62 


24 


/ 


/ 


+19.14 


10.12 


+3.73 


25 


/ 


/ 


/ 


-7.17 


-5.55 


26 


/ 


/ 


/ 


/ 


+3.26 


34 


/ 


/ 


-179.31 


-172.61 


-159.03 


35 


/ 


/ 


/ 


+11.78 


+3.26 


36 


/ 


/ 


/ 


/ 


-5.05 


45 


/ 


/ 


/ 


-187.04 


-180.62 


46 


/ 


/ 


/ 


/ 


+11.58 


56 


/ 


/ 


/ 


/ 


-191.51 



In Table II, we provide details on the rotation angles 
and the magnetic moments of the non-collinear states of 
the even-numbered chains. In addition, we calculated the 
energy differences AE between the non-collinear state 
and the ferrimagnetic solution for even- numbered chains, 
showing that the non-collinear states are energetically 
favored. 



TABLE II. Ab initio results for the nanochains: size and 
angle of the magnetic moments as well as total energy dif- 
ferences between the non-collinear and antiferromagnetic 
solutions. In every pair of wire atoms connected by a "-" sign, 
the azimuthal angles (j) are equal to O'^ and 180^ respectively, 
while the magnetic moments and rotation angles 6 are the 
same. 



Length 
(adatoms) 



Adatom 6>(^) M (/is) 



(meV/adatom) 



1-2 



77 



3.53 



1-4,2-3 67, 85 3.52, 3.29 



-2.26 



-5.59 



1-6, 2-5 57, 94, 3.51, 3.30 -6.33 
3-4 76 3.27 



3 



Zero-point energy of a tetramer 

The total magnetic energy of a tetramer considering 
only next neighbor exchange interactions and uniaxial 
anisotropy is given by: 

^4(^1,^2,01,02) 

= -2 Ji2 cos((9i + O2) - J23 cos(2(92) 
- 2 J/i cos(6>i) - 2 J/2 cos(i92) 
+ 2i^(sin2(^i) cos2((/)2) + sin\02) cos^{h))- 



(1) 



where J12 and J23 are the exchange constants given in 
Table 1. J/i and J/2 account for the ferromagnetic ex- 
change between the substrate and the inner and outer 
Mn atoms, respectively. From the DFT calculations we 
obtain the values J/i =90 meV and J/2 = 94 meV. The 
magnetic anisotropy energy K per atom is taken from 
the value of the dimer. The orientation of the magnetic 
moment of the inner atoms is described by ^1, 01 and 
by ^2, 02 for the outer atoms. The angles in the ground 
state minimize the energy given in equation 1 and are 
found to be: 

^10 = 65° 620= 88° 
010 = 0° 020 = 180° 

where the additional subscript stands for ground state. 
Note that this approach leads to results similar to our 
DFT calculations. The zero-point energy of the tetramer 
Eo4 is found to be [12]: 



Eo4 = 
with 
b = 



-huj 




(3) 



^6>i6>i^(/)i(/)i 



M|sin2(^2o) 

ESnSn — EnL 



Es,6, Es 



where g ^ 2 \^ the g-factor and fiB the Bohr magneton. 
Efy^j^ is the second derivative of ^4 with respect to a and 
P evaluated at the optimum angles, i.e.. 



d?E 



dced/3 



As ^^202 = ^^101 = for 01 = 010 and 02 = 02o, this 
derivative does not explicitly appear in (3). Mi and M2 
are the magnetic moments obtained from the DFT cal- 
culations. The lowest possible zero-point energy is thus 
given by 



Eo4 = 5.2 meV, 



(4) 



which is again larger than the barrier for the rotation 
in the 0-direction, given by \2K\{sm^{0i) + sin^(6>2)) = 
1.1 meV. 
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